function base_val = nature_base_func(X,K)
%X=reshape(X,[length(X),1]);
base_val=zeros(length(X),K+1);
for iiK=0:K
base_val(:,iiK+1)=X.^iiK;
end
end
